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Abstract 

The meaning of time in an open quantum system is considered under the assumption that 
both, system and environment, are quantum mechanical objects. The Hamilton operator of the 
system is non-Hermitian. Its imaginary part is the time operator. As a rule, time and energy 
vary continuously when controlled by a parameter. At high level density, where many states 
avoid crossing, a dynamical phase transition takes place in the system under the influence of the 
environment. It causes a dynamical stabilization of the system what can be seen in many different 
experimental data. Due to this effect, time is bounded from below: the decay widths (inverse 
proportional to the lifetimes of the states) do not increase limitless. The dynamical stabilization 
is an irreversible process. 
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I. INTRODUCTION 



The problem of time is considered in very many papers starting from the very beginning 
of quantum mechanics. For example, the derivation of the uncertainty relation between 
space and momentum can be found in every textbook on quantum mechanics. However, the 
uncertainty relation between energy and time could not be derived convincingly up to now. 
According to Pauli the reason is that energy is bounded from below while this is not 
the case for time and, furthermore, time varies continuously in contrast to energy. A critical 
consideration of the discussions of the time energy uncertainty relation is given in [2]. The 
authors point out that energy E characterizes the quantum system while time t is measured 
by an external clock. 

n 

This point of view is followed up in recent studies [3]: it is important in the application 
of a given time energy uncertainty relation to state precisely what kind of measurement is 
being made and to specify accordingly the meaning of the At involved, e.g. does it refer to 
the accuracy of measurement, to the duration of measurement or perhaps to the lifetime of a 
decaying state. Also the derivation of the time dependent Schrodinger equation is considered 
4j: This derivation, unlike those presented in quantum mechanics textbooks but in the spirit 
of Schrddinger's original approach to the problem, acknowledges that time enters quantum 
mechanics only when an external force on the quantum system is considered classically. 
Starting from a fully time-independent formulation of quantum mechanics, it is possible 
systematically to derive the time- dependent Schrodinger equation for a quantum system in 
the approximation that the environment is treated semiclassically. The derivation relies on 
the assumption of a closed object comprising system plus environment 4]. 

In all these considerations the Hamiltonian of the system is assumed to be Hermitian, 
and the lifetime of a decaying state is calculated perturbatively. The situation is another one 
when both, system and environment, are quantum mechanical objects, i.e. the environment 
is the continuum of scattering wavefunctions into which the system is embedded. In such a 
case, the Hamilton operator of the system is non-Hermitian and the lifetime of a decaying 
state is calculated non-perturbatively |5j. 

It is the aim of the present paper to consider the meaning of time in an open quantum 
system described by a non-Hermitian operator the eigenvalues of which are complex. The 
basic equations are taken from the review |5[ . It turns out that a natural definition of time 
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in an open quantum system is related to the decay widths of the states and, furthermore, to 
the dynamical stabilization occurring in the system under the influence of the environment. 
Time, defined in such a manner, is bounded from below. 

The paper is organized in the following manner. In section |TT] it is shown why in an open 
quantum system with many levels the non-Hermitian operator H e g appears and how it looks 
like. Using this model, the properties of discrete and narrow resonance states are sketched 
in section IHIl In both cases, exceptional points play an important role. Since there is much 
confusion in the literature on these singular points, their properties are sketched in the 
appendixes \M to [E] from a unified point of view and compared to the results of experimental 
studies. The relation of the imaginary part of the eigenvalues of H c $ to time is considered 
in section HVl while in the following section |V] the decay rate is calculated. The phenomenon 
of a dynamical phase transition in an open quantum system is explained in section I VI I In 
section IVHt a few examples of dynamical phase transitions, that are found experimentally, 
are given. In any case, they cause a dynamical stabilization of the system. Based on these 
results, it is possible to define time in an open quantum system and to discuss its properties 
(section IVIHI) . The results are summarized in the last section HXl 



II. NON-HERMITIAN HAMILTON OPERATOR OF AN OPEN QUANTUM 
SYSTEM 

The definition of an open quantum system used in the present paper is the following: a 
quantum system is considered to be open when it is embedded into the environment of the 
continuum of scattering wavef unctions. According to this definition, the system is localized 
while the environment is extended infinitely. This type of embedding of a quantum system 



into an environment always exists. The environment may 



by a laser in the case of atoms (examples are shown in {(5]). It can however not be deleted 
completely. When allowed by the energy of the system (and no special selection rules hold), 
most states of the system decay into the continuum of scattering wavefunctions and have a 
finite lifetime. Otherwise they are discrete. 

A method being very suitable for the description of this situation, is a projection operator 
formalism with projection onto system and environment, respectively. In the following, 
some basic equations of such a formalism are given, for details see the review {J. First the 



)e changed by external forces, e.g. 
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Schrodinger equations in the two subspaces (with the Hamiltonian H p and H c , respectively) 
have to be solved and the projection operators Q and P have to be defined, 

N 

(H B -Ei)^ = o -+ Q = £|$fx<&f| (i) 

k=l 
C poo 

(H c -E)^ = -> P = W dE\g)(g\. (2) 

c=l J -°° 

The $f are the wavefunctions describing the N discrete states of the closed many-particle 
system, while the £,f are the scattering wavefunctions of the environment consisting of C 
continua (£,f is written instead of £ ( f( + ' ) f° r convenience). The Schrodinger equation in the 
whole function space with discrete and scattering states reads 

(H fun - E) tff = (3) 

with the Hermitian operator # MI = H QQ + H QP + H PQ + H PP and Q + P = 1 where 
Hqq = QHQ and so on. The coupling matrix elements between system and environment 
are 

ll = ^{^\H QP \^). (4) 
The solution of the full problem (jHJ) is 

«f = g +£n>- < * t ' g< ^ ltf> (5) 

fc=l ^ - ^fc 

with the non-Hermitian operator 

-f^cff = Hqq + Hq P G p H P q = H B + VbcG^c Vcb (6) 

and 

(H eS -z k )$ k = 0; z k = E k - l -Y k (7) 

after diagonalization. Here E k = E k (E) and T k = T k (E) are, respectively, the position in 
energy and the decay width (inverse lifetime) of the state k at the energy E, see equations 
(ED and ©. 

Position and width of the resonance state k are energy independent numbers E k = El 
and T k = T p k only when the state k is not overlapped by another resonance state and, 
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furthermore, it is far from a decay thresho 
usually, from the poles of the S matrix, see 



d. These numbers can be obtained also, as 
7|. The functions Zk = Zk(E) describe, in any 



case, the line shape of resonances correctly. Further, Gp~^ = P(E — Hpp)~ l P is the Green 
function in the P subspace and Vt k = (l + G^ Hpq) is the wavefunction of the resonance 
state. The eigenf unctions of H e{i are biorthogonal, (<&£|<&j) = S^i and ($ fc |$ fc ) = A k > 1. 

The Hamiltonian H c q consists formally of a first-order and a second-order interaction 
term. The second-order term via the continuum determines the dynamics of the system at 
high level density. It leads to the principal value integral 

1 C K 7° 7° 

c=l 



Re mH^-EfS^ = jTdtf ^| (8) 

and the residuum 



Im <*f \H eS \$f > = -\ f>°7? c . (9) 

Z c=l 

When i = j, the expression Re($f \H e e\^f ) gives the shift in energy of the state i due 
to the interaction between the state and the large system it is part of, i.e. due to the 
interaction of the state % with the environment (subspace P). This is the self-energy of the 
state which is analog to the Lamb shift known in atomic physics. When i ^ j, OH]) describes 
the energy shift of the state i due to its coupling to another state j ^ i via the environment 
(continuum of scattering wavefunctions) . These couplings to the different states j cause, 
at high level density, collectively contributions to the energy shift of the state i. In atomic 
physics, contributions of such a type are studied newly experimentally and called collective 
Lamb shift 0, J 

The energy window coupled directly to the continuum is e c < E < e' c . In nuclei e' c — > oo 
while e c denotes the lowest threshold for emission of a particle. The energy shift = 
Ek — Efr, including the corresponding corrections arising from the coupling of different states 



via the continuum, can not be simulated by two-body forces [10|. It is a global property, see 



11] , and is not contained in any standard calculation with a Hermitian Hamilton operator 
for the many-body problem. 

The method sketched here for the description of an open quantum system is not the 
only one. The advantage of the explicit consideration of H e g in the model described above, 
consist, above all, in the fact that the many-body problem in the Q subspace has to be solved 
only once, since it is energy independent, see (JT|). At low level density, H cS w H B = Hqq, 



the second-order term can be treated perturbatively, i.e. the Hamiltonian H e g is almost 
Hermitian, corresponding to the assumption of standard quantum physics. At high level 
density, however, the second-order term can not be treated perturbatively. It induces a 
global mixing of the states according to (jSJ) [111 ], and deviations from standard quantum 
physics will occur. They appear most clearly when the number of decay channels (continua) 
is small, especially for C = 1. 



III. DISCRETE AND NARROW RESONANCE STATES 

The eigenvalues of the effective Hamiltonian H e g, equation fl2]), may be real or complex. In 
the first case, the eigenstates are discrete states while they are resonance states in the second 
case. The boundary conditions are different for the two different types of states. In both 
cases, the trajectories of the eigenvalues traced as a function of any control parameter avoid 
crossing, usually. The corresponding crossing points (called mostly exceptional points) of two 
eigenvalue trajectories can be found by analytical continuation. They play an important role 
for the dynamics of open quantum systems. The eigenvalues and eigenfunctions at and in 
the neighborhood of an exceptional point (including the experimental proof of their main 
features) are given in appendix [A] and [Bj respectively. 

Resonance states are coupled directly to the continuum, e c < E k < e' c . The eigenvalues 
Zk of H e Q are complex, generally. According to ([5]), the scattering wavefunction int inside 
the system can be represented in the set of the biorthogonal eigenfunctions of if e ff, i-e. 

*f - - £ <*ijf« r 5> (10) 

fc=l z k 

The expression (<&j£|<&j) is a complex number, and a consistent normalization of the $fc 

requires Im($£|$fc) = 0. This corresponds to some rotation such that the phases of the 

eigenfunctions relative to one another are not rigid (when traced as a function of a certain 

fl 

parameter). Instead, it holds < p < 1 (for details see |5|) for the phase rigidity 



where 9 is the rotation angle (see appendix O for the definition of the phase rigidity r k 
in the case of a two- level system). Only at low level density and far from avoided level 
crossings, ~ and p « 1 such that standard Hermitian quantum physics is 



a good approximation. Otherwise, spectroscopic redistribution processes take place, p < 1, 
the Schodinger equation has a nonlinear source term (see appendix |D]) and time reversal 
symmetry is broken (see appendix [E]). Here, the description of the system by standard 
Hermitian quantum mechanics breaks down. 

The breakdown of standard quantum mechanics can be understood in the following man- 
ner. When 1 > p > the states with the wavefunctions int avoid crossing (when controlled 
by a parameter) and become mixed globally such that a few states of the system can align 
hierarchically (i.e. step by step with the scattering states ^ of the environment. As 
a result, their decay widths become large. Full alignment is reached for p = 0. The 
alignment of a few states tyf int with the channel wavefunction £^ occurs by trapping other 
resonance states, i.e. by (partial or complete) decoupling them from the environment. This 
is nothing but width bifurcation: the widths of a few states become large while the widths 
of the other ones become small by varying the control parameter. This scenario occurs in 
the vicinity of exceptional points j^l, Q]. When p < 1, the system can not be described 
perturbatively, and standard Hermitian quantum physics fails (compare Appendixes [C 
for the two-level system). Similar results are obtained by using other methods, e.g. [13 

The influence of the continuum of scattering wavefunctions (P subspace) onto the discrete 
states of the system with energies Ek beyond the energy window e c < E < e' c , seems to be 
much less important. For discrete states H e g is non-Hermitian, but the eigenvalues Zk = E^ 
are real, the eigenfunctions are orthogonal, = dik, and the phase rigidity is p = 1. 

However, Ek ^ E^ . According to (IE]), the energy shift = Ek — Ej? ^ is caused by 
the coupling of the state k to the environment, i.e. by the embedding of the system into 
the continuum of scattering wavefunctions. By this, many-body forces are induced in the 



system 



101 ] . Discrete states avoid crossing and, at a critical value of the control parameter, 



the two states are exchanged as known for about 80 years 17J. The only difference to the 
avoided crossing of resonance states is that discrete states never cross. The corresponding 
crossing point can be found only by analytical continuation into the continuum [18j ]. 

When the level density is high, many discrete levels avoid crossing (when traced as a 
function of a control parameter). For illustration let us consider an A particle system. In 
this case, the induced many-body forces cause a global mixing of the discrete states in a 
finite parameter range, see Finally, an aligned discrete state is formed with the structure 
bound particle + (A-l) particle residual system which corresponds to the structure of the 
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decay channel unbound particle + (A-l) particle residual system (the quantum numbers of 
particle and residual system are the same in both cases and (A-l) denotes the number of 
particles of the localized residual system after emission of one particle into the continuum). 
This discrete state is the analog to an aligned resonance state above particle decay threshold. 
The only difference between these two states is that the energy of the preformed aligned 
discrete state is too small and does not allow the emission of one particle from the system, 
while the decay width of the aligned resonance state is large, corresponding to a short lifetime 
of this state. 

It should be mentioned that energy conservation is not the only source for an eigenvalue of 
H eS to be real. Im(z fc ) = — T k /2 = is possible also due to selection rules (according to the 
corresponding quantum numbers of the states), or because of width bifurcation appearing 
at high level density (causing the so-called bound states in the continuum [5]). Another 
source for real eigenvalues of a non-Hermitian operator app ears in PT symmetric systems 



(where P and T denote parity and time, respectively) [19] and can be observed in optics 
due to the formal equivalence of the optical wave equation in PT symmetric optical lattices 
to the quantum mechanical Schrodinger equation 20]. The relation of these results to 



the properties of open quantum systems as discussed above, is considered in 
appendix [A] 



12], see also 



IV. THE IMAGINARY PART OF THE EIGENVALUES OF THE OPERATOR H eS 

In order to study the physical meaning of Im(zk) = — r^/2, equation (jTJ), the behavior of 
three neighboring resonances in a two-dimensional quantum billiard connected to a single 



2l] . A measurable quantity derived from the 



waveguide is investigated theoretically in 
reflection coefficient R{E) is the Wigner Smith time delay function 

It is the time the wave spends inside the billiard. The energies E% and widths Y\ of the 
resonance states can be found from the poles of the function R(E) analytically continued 



into the 
In 



ower complex plane. 



21] , contour and surface plot of ln{r w ) and the motion of the corresponding three 



resonance poles by varying the coupling between the waveguide and the resonator are cal- 



S 



culated. At weak coupling to the waveguide, the three resonances are seen clearly. As the 
coupling to the waveguide increases, the lifetimes of all three states decrease, as expected. 
As the resonances start to influence one another, the states attract each other in energy, 
two of them become trapped while the third one becomes short-lived. At further increasing 
coupling, the lifetimes of the two trapped resonance states increase, contrary to expectation. 
The lifetime of the short-lived state is, at large opening, so short that it practically disap- 
pears when plotting the time delay. The motion of the poles is reflected in the time delay 
function. The calculations have shown further that the interference between the resonance 
states leads to a mixing of their wavefunctions with respect to the eigenfunctions of the 
closed resonator (defined by decoupling the resonator from the waveguide). 

Some years ago, the dynamics of resonance states is studied experimentally by means of 
a flat microwave resonator connected to a waveguide where the coupling strength between 
resonator and waveguide can be varied by hand j^ . In this experiment, the microwaves 
enter the billiard through a slit, the opening of which can be varied. The motion of the 
resonance poles as a function of the opening of the slit is traced starting close to the real 
axis and following them into the region of overlapping resonances. The results verify the 
resonance trapping effect discussed above. This experimental proof does not depend on any 
model assumptions. Meanwhile, the resonance trapping effect has been investigated and 
verified in many other theoretical and experimental studies performed on different systems 



(see section IVHl and review {5]). 

As a result of these studies, the physical meaning of the imaginary part Im(zk) of the 
eigenvalues of the non-Hermitian Hamilton operator H c g is directly related to the time r w . 
This relation holds also at high level density where the system as a whole is dynamically 
stabilized (see section IVTl). 



V. DECAY RATE AT HIGH LEVEL DENSITY 

The time dependent Schrodinger equation for the scattering wavefunctions tyf int inside 
the system reads 

H M1 *f int (i) = ih^ Ut) (13) 



9 



with H inl1 defined in and int in ({TO]) . The right solutions |\I/f int ) may be represented 
by an ensemble of resonance states k that describes the decay of the localized part of the 
system at the energy E, 

N 

l*f !*(*)> = e- iHeBt/% l*?int(*b)> = E e^^cf fe |$ fc ) (14) 

fc=i 

with cf fc = (^|F|Cf)/(-B - z k ) according to (TTUT) . By means of ([14"]) and the correspond- 
ing expression for the left solution of f[T3"j) . the population probability (\I/ c (t)|\I/ c (t)) = 
Ea c c 2 fe e" rfc * /?i with the energy averaged values c ck can be defined. The decay rate reads [23] 

8 1 V, r, r 2 p- T kt/h 

U*) = ~ ¥t ^ <*? int (*)l*f int (*)> = j- \ c £ k e -v kt/ n ■ (15) 

The decay properties of the resonance states can be studied best when their excitation takes 
place in a time interval that is very short as compared to the lifetime t\ of the resonance 
states. In such a case, no perturbation of the decay process by the still continuing excitation 
process will take place. 

For an isolated resonance state k, f|15j) passes into the standard expression 
fcgr(t) — > k k = Tk/h. In this case, the quantity k^ is constant in time and corresponds 
to the standard relation r k = h/Tk with = It describes the idealized case with 

exponential decay law and a Breit-Wigner resonance in the cross section. 

Equation 03 descnbes, however, the decay rate also in the regime of overlaps res- 
onances [23| . The overlapping and mutual influence of resonance states is maximal at the 
avoided (and true) crossing points in the complex plane where two eigenvalues and zy of 
the effective Hamilton operator H e g coalesce (or almost coalesce). Nevertheless, the decay 
rate is everywhere smooth as can be seen also directly from ( 1T51) . This result coincides 
with the general statement according to which all observable quantities behave smoothly at 
singular points. 

An interesting result is the saturation of the average decay rate k av in the regime of 
strongly overlapping resonances. According to the bottle-neck picture of the transition state 



theory, it starts at a certain critical value of bound- continuum coupling 2% l25|. As has 
been shown in |26| |. this is caused by width bifurcation since the definition of an average 
lifetime of the resonance states is meaningful only for either the long-lived states or the 
short-lived ones. The widths Y k of the long-lived (trapped) states are almost the same for 
all the different states k, i.e. T av w T k for all long-lived resonance states |5|. It follows 

10 



therefore k av ~ T av /h from (fl5l) . According to the average width T av , the average lifetime 
of the long-lived states can be defined by r av = l/k av . Then r av = h/T av . That means, 
the basic relation between lifetimes and decay widths of resonance states holds not only for 
isolated resonance states but also for the narrow (trapped) resonance states, i.e. for T av and 

7~ecv ■ 



VI. RESONANCE TRAPPING AND DYNAMICAL PHASE TRANSITIONS 



Some years ago, the question has been studied 11] whether or not the resonance trapping 
phenomenon is related to some type of phase transition. The study is performed by using 
the toy model 

Hl° y = H + iaVV + (16) 

where H° and VV + are Hermitian, V is the coupling vector of the system to the environment 
and the parameter a simulates the coupling strength between system and environment. The 
calculations are carried out for the one-channel case and with the assumption that (almost) 
all crossing (exceptional) points accumulate in one point [271 ] . The control parameter a is a 
real number. It has been found that resonance trapping may be understood, in this case, as 
a second-order phase transition. The calculations are performed for a linear chain consisting 
of a finite number N of states. The state in the center of the spectrum traps the other 
ones and becomes a collective state in a global sense: it contains components of almost all 
basic states of the system, also of those which are not overlapped by it. The normalized 
width T /N of this state can be considered as the order parameter: it increases linearly as a 
function of a, and the first derivative of Tq/N jumps at the critical value a = a cv . The two 
phases of the system differ by the number of localized states. In the case considered, this 
number is N at a < a cr , and N — 1 at a > a cr . 

Much more interesting is the realistic case with the Hamiltonian ([6]). In this case, trapping 
of resonance states occurs in the regime of overlapping resonances hierarchically, i.e. one 
by one {^. The crossing points do not accumulate in one point, but are distributed over a 
certain range of the parameter: a dynamical phase transition takes place in a finite parameter 
range inside the regime of overlapping resonances. It can therefore be observed [5|. Also in 
this case, almost all resonance states are involved in the phase transition of the system and, 
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furthermore, the number N of localized states is reduced. 

The dynamical phase transition taking place in the system at high level density causes 
finally a dynamical stabilization of the system: the system consisting of only the localized 
long-lived states beyond the phase transition is more stable than the system below the 
dynamical phase transition in spite of the stronger coupling between system and environment 
beyond the phase transition. The reason is the following: first the widths bifurcate at 
high level density and then the state with the shortest lifetime is ejected. The dynamical 
stabilization is a global effect to which all states contribute collectively by aligning one of 
the states (step by step) with a decay channel. In this manner, the sum of the decay widths 



of the states of the system is reduced [28[, and the system is stabilized. 



VII. DYNAMICAL STABILIZATION OF DIFFERENT QUANTUM SYSTEMS 
IN EXPERIMENTAL DATA 

A. Phase lapses 

More than 10 years ago, in experiments on Aharonov-Bohm rings containing a quantum 
dot in one arm, both the phase and the magnitude of the transmission amplitude T = \T\ e l/3 



of the dot are extracted [29|. The results obtained caused much discussion since they do 
not fit into the standard understanding of the transmission process. As a function of the 
plunger gate voltage V g , a series of well-separated transmission peaks of rather similar width 
and height has been observed in many-electron dots and, according to expectations, the 
transmission phases (3(V g ) increase continuously by ir across every resonance. In contrast to 
expectations, however, always jumps sharply downwards by ir in each valley between any 
two successive peaks. These jumps called phase lapses, were observed in a large succession 
of valleys for every many-electron dot studied. Only in few-electron dots, the expected so- 
called mesoscopic behavior is observed, i.e. the phases are sensitive to details of the dot 
configuration. The problem is considered theoretically in many papers over many years 



without solving it convincingly, e.g. 



30). 



In [3JJ, the generic features of phase lapses in the inelastic cross section are studied 
by using the toy model ( fl6i) for the non-Hermitian Hamilton operator. According to the 
results of these calculations, the universal features observed in the phase lapses at high level 
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density, in contrast to the mesoscopic features at low level density, may be considered to 
be a hint at a dynamical phase transition. The transition occurs by controlling the system 
from low to high level density simulated in the calculations by means of a. In accordance 
to this picture, only the resonance states at low level density show individual spectroscopic 
features. At high level density, the observed resonances arise from trapped states. They show 
level repulsion, have vanishing spectroscopic relation to the open decay channels (i.e. small 
decay widths), and phase lapses appear. It follows further, that any theoretical study on 
the basis of conventional Hermitian quantum physics is unable to explain the experimental 
results convincingly. More accurate calculations on the basis of ([7]) for the non-Hermitian 
Hamiltonian H e g are performed recently and compared with the experimental data. 



B. Spin swapping operation 

A swapping gate in a two-spin system exchanges the degenerate states | f, \) and | |, t)- 
Experimentally, this is achieved by turning on and off the spin-spin interaction b that splits 
the energy levels and induces an oscillation with a natural frequency u>. An interaction Ti/tse 
with an environment of neighboring spins degrades this oscillation within a decoherence time 
scale t^. The experimental frequency u> is expected to be roughly proportional to b/h and the 



decoherence time proportional to t$e- In 13j . experimental data are presented that show 



drastic deviations in both u and from this expectation. Beyond a critical interaction 
with the environment, the swapping freezes and the decoherence rate drops as l/r^ oc 
(b/h) 2 TsE- That means, the relaxation decreases when the coupling to the environment 
increases. The transition between these two quantum dynamical phases occurs when u oc 



(b/h) 2 — (k/rsE) 2 becomes imaginary (where k depends only on the anisotropy of the 
system-environment interaction, < k < 1). The experimental results are interpreted by 
the authors as an environment ally induced quantum dynamical phase transition occurring 
in the spin swapping operation |l3| . 

Further theoretical studies within the Keldysh formalism showed that is a non-trivial 
function of the system-environment interaction rate t S e, indeed: it is oc 1/t S e at low 
tse (according to the Fermi golden rule) but l/r^ oc tse at large t<?b. This theoretical result 
is in (qualitative) agreement with the experimental results. In [13J, the dynamical phase 
transition in the spin swapping operation is related to the existence of an exceptional point. 
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The dynamical phase transition observed experimentally in the spin swapping operation 
and described theoretically within the Keldysh formalism shows qualitatively the same fea- 
tures as the dynamical phase transitions discussed in the present paper on the basis of the 
resonance trapping phenomenon (width bifurcation). 



C. Loss induced optical transparency in complex optical potentials 

The optical wave equation for complex PT symmetric potentials is formally equivalent 



to the quantum mechanical Schrodinger equation [33[. One expects therefore that PT sym- 
metric optical lattices show a behavior which is qualitatively similar to that discussed for 
open quantum systems. 

Experimental studies showed, indeed, a phase transition that leads to a loss induced 



optical transparency in specially designed non-Hermitian guiding potentials 20]: the output 
transmission first decreases, attains a minimum and then increases with increasing loss. 
The phase transition is related, in these papers, to PT symmetry breaking. In a following 



theoretical paper [34|, the Floquet-Bloch modes are investigated in PT symmetric complex 
periodic potentials. As a result, the modes are skewed (nonorthogonal) and nonreciprocal. 
That means, they show the same features as modes of an open quantum system under the 
influence of exceptional points. A detailed discussion of this analogy is given in [l2j. 



D. Dicke superradiance and subradiance in optics 



The assumption that the probability of a given molecule to emit a photon may be consid- 
ered to be independent of the states of the other molecules is justified only when the distance 
between the molecules is large. Generally, all the molecules are interacting with the common 
radiation field and the spontaneous emission takes place coherently. Dicke 35[ was the first 
who considered the coherence in spontaneous radiation processes and, as a consequence, the 
formation of the so-called superradiant state. The collective coupling of the atoms via the 
radiation field leads also to a substantial radiative shift of the transition energy, the so-called 
collective Lamb shift. This effect is recently proven experimentally {§]. It will allow us to 
probe aspects of quantum electrodynamics in relatively low-energy experiments js]. 

By using the simple non-Hermitian Hamilton operator f|T6|) . the formation of a superradi- 
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ant state can be controlled by the coupling strength a between system and environment 36 ]. 
At large a, the short-lived superradiant state is formed together with long-lived subradiant 
(trapped) states. The corresponding phase transition, called superradiance transition by 
the authors, is nothing but the dynamical phase transition discussed in section IVII In both 
cases, the same Hamiltonian with frozen internal degrees of freedom is used. In section IVII 
the relation to exceptional points is considered. 

Using an effective non-Hermitian Hamilton operator describing photon-mediated atomic 
dipolar interactions, calculations for the Dicke superradiance in atomic gases are performed 



recently 



37J. The calculations show all the features characteristic of a dynamical phase 



transition. Only the notations used are different: for example 'disorder' used in [37] cor- 
responds to 'individual spectroscopic properties of the states' in the many-body problem 
(sections UTI and [TTT]) . Also in these calculations, the crossover from the uncorrelated to the 
correlated behavior appears at a critical value of a certain parameter. Interesting is the 
scaling behavior of the escape rates of photons propagating in a 3D atomic gas when traced 
by means of a certain parameter. For small parameter values, cooperative effects are negli- 
gible, and photons are emitted in spontaneous and incoherent processes. For larger values 
of this parameter, cooperative effects set in and the Dicke superradiance appears. The cor- 
responding subradiance is considered in 



It extends the lifetime of the excitation to 



many times the natural lifetime of a single atom and is, therefore, doubtless interesting in 
quantum information science. 



VIII. TIME IN AN OPEN QUANTUM SYSTEM 

The results of sections IIVI to IVHI show that resonance states with very large widths, 
corresponding to very short lifetimes, do not exist. Due to the avoided crossings of the 
eigenvalue trajectories and the accompanying reduced phase rigidity p < 1, most resonance 
states are dynamically localized at high level density and have a finite non-vanishing lifetime. 
The following conclusions can be drawn. 

(i) Energy and time are determined by the eigenvalues of one and the same operator, namely 
of H e ff. The time operator is the non-Hermitian part of if e fr while the energy operator 
corresponds, as well known, to the Hermitian part of H e g. 

(ii) Not only time varies continuously (as well known) but also energy does so, since most 
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states of an open quantum system are resonance states with a non-vanishing width. 

(iii) Time is bounded from below in a similar manner as energy. 

(iv) Discrete states correspond to t — > oo. This relation is analog to the assumptions of 
standard quantum mechanics described by Hermitian operators. 

The two operators Ke(H e g) and Im(if e fj) do not commute when the number C of channels 
is different from the number N of states, as can be seen easily from ( fl6|) . Usually C <C N. 
According to the points (i) to (iv), the main argument by Pauli [l| against the derivation of 
the time energy uncertainty relation does not exist when an open quantum system with the 
non-Hermitian Hamilton operator H e $ is considered. 

The concept of time considered in the present paper is characteristic of an open quantum 
system in the same manner as energy. It has a physical meaning only for £& < t < oo where 
tb is the time at which width bifurcation creates the dynamic stabilization of the system. 
The time t defined in this manner is a measurable quantity. Note t is more sensitive to 



parameter variations than E (see 



39( and 53J for numerical examples). 



The dynamical stabilization of the system at high level density appears as a counter- 
intuitive process when time is considered as a parameter that may be varied continuously 
between — oo and +00 and is not characteristic of the system. In contrast to this, the finite 
value tb relies on the fact that time is a value characteristic of the system. The dynamical 
stabilization occurs in consequence of the fact that the resonances will never really overlap. 
The only way to achieve this is to accumulate almost all coupling strength between system 
and environment onto one state (in the one channel case) while the remaining states become 
stabilized (long-lived) 40(]. The short-lived state created in this manner, is aligned to the 
continuum of scattering states and does not belong to the set of localized states. Beyond tb, 
the system differs from the original one: the number of states is reduced and their individual 
spectroscopic features are lost. Instead, cooperative effects are important (which may be 
disturbed only by coupling the system to another channel). Mathematically, the redistri- 
bution rests on the existence of exceptional points and the related phenomenon of avoided 
level crossings. It occurs under the influence of the environment into which the system is 
embedded. 

Experimentally it is possible to form spatially remote discrete states on separate quantum 



point contacts and to allow an interaction between them via a common continuum 



42| . The 



results show that the continuum supports an effective interaction between the two states 
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which is mediated by the continuum and is highly robust: the detector exhibits two Fano 
resonances due to the two different bound states. These studies, revealing a pronounced 
avoided level crossing, show clearly that the continuum affects the microscopic structure of 
bound states in the two quantum point contacts. 

In the examples discussed in section IVIH the time i& appears in a natural manner. It 
causes measurable effects. For example, the transmission through a small quantum system 
is enhanced in the parameter range in which the redistribution in the system takes place 



4=11 ] . Interesting are the very stable whispering gallery modes in a small quantum system. 



They are partly aligned to the scattering wavefunctions and their lifetimes are shorter than 
those of the other states. They cause an enhancement of the transmission, and the system 
becomes almost transparent (p — > 0) under these conditions. 

The time tb determines also the brachistochrone problem 5| which consists of finding 
the minimal time for the transition from a given initial state to a given final state of the 
considered system. At high level density, the individual resonance states can no longer be 
identified. Here p < 1, and the wavefunctions of some states of the system are partly aligned 
to the scattering wavefunctions such that the time for the transition from a given initial state 
to a given final state may be radically shortened. However, the time for traveling through 
the system does never vanish. It is bounded from below since it can not be smaller than the 
time corresponding to the transparency of the system 5|. 

The non-adiabatic processes found recently when cycling exceptional points 43], are 
surely related to the finite time tb below which time loses its meaning in the system consid- 
ered. The cycling crosses regions with fundamental different time concepts. 



IX. SUMMARY 



In the present paper, the meaning of time in an open quantum system is considered. 
System as well as environment are quantum mechanical objects with the consequence that 
the Hamilton operator H c q of the system is non- Hermit ian. The real parts of the eigenvalues 
of H e s provide the energies of the states while the imaginary parts of them are related to 
their lifetimes. Time t in the open quantum system is defined by the lifetimes of the decaying 
states. 

In an open quantum system described by a non-Hermitian operator, the main objection 
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[l| to the derivation of the time energy uncertainty relation does not occur. The time 
operator appears in a natural manner together with the energy operator. Energy and time 
are bounded from below and, furthermore, vary continuously as a function of a parameter, 
as a rule. 

The open quantum system is reversible at low level density where the levels are far from 
one another and the system can be described, to a good approximation, by a Hermitian 
operator. At high level density, however, the levels avoid crossing and irreversible processes 
determine the redistribution processes taking place in the system. These irreversible pro- 
cesses are caused by nonlinearities in the Schrodinger equation of the open system in the 
vicinity of avoided level crossings (and exceptional points, respectively) due to the coupling 
of the states via the environment. A dynamical stabilization of the system occurs, at the 
time tb, under the influence of the environment, tb is the lowest value of t. 

Further experimental as well as theoretical studies are necessary. Above all, the time 
energy uncertainty relation has to be derived in a mathematical convincing manner for the 
general case with discrete and narrow resonance states of a many-particle system. Experi- 
mentally, the influence of a third state onto the mixing of two states in the neighborhood 
of an avoided level crossing should be studied. As discussed in appendix [El this is the basic 
process of the dynamical stabilization taking place in the system at high level density. 

Appendix 

Appendix A: The eigenvalues of a non-Hermitian 2x2 operator 

Let us consider the Hamiltonian 



with the energies q (i — 1, 2) of the two states and the interaction u between them. The e; 
are assumed to contain the corrections due to the coupling of the state i to the environment. 
The eigenvalues are 



The levels repel each other in energy according to the value Ke(Z) while the widths bifurcate 
corresponding to lm(Z). The two eigenvalue trajectories cross when Z = 0, i.e. when (ei — 




(Al) 




(A2) 
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t2)/2uj = ±i. At the crossing points, called mostly exceptional points, the two eigenvalues 
coalesce, E\ = e<i = eq. In the vicinity of the crossing points, the dependence of the 
eigenvalue trajectories on the parameter is more complicated than far from them: the two 
levels approach each other in energy and the widths become equal so that Re(ei) -H- Re^) 
and Im(ei) -H- Im(e 2 ) at the crossing point. 

When if is a Hermitian operator, the unperturbed energies q of the states and the 
interaction u between them are real. According to (IA2I) . the two (real) eigenvalue trajectories 
Ei(a) = ej(a) cannot cross (for u^O) when traced as a function of a certain parameter a. 
Instead, they avoid crossing. The Active crossing point is called usually diabolic point. The 
topological structure of this point is characterized by the Berry phase [44j which is studied 
theoretically and experimentally in many papers. 

The situation is another one when if is a non-Hermitian operator. In such a case, the 
unperturbed energies e, and also the interaction oj are complex, usually. The states can 
decay, in general, and the two eigenvalues ( 1A2I) can be written as 

£1,2 = ei, 2 - - 7i,a (with 71,2 > 0) . (A3) 

The widths 7$ are proportional to the inverse lifetimes r^ 1 of the states, i — 1,2. The 
two eigenvalue trajectories £i(a) may cross according to (1A2I) . and the crossing point is 
an exceptional point in agreement with the definition given in 45j. The topological phase 
of the exceptional point is twice the Berry phase 5|. This theoretical result is proven 
experimentally by means of a microwave cavity 49 ]. According to (1A2|) . Re(Z) causes 
repulsion of the levels in energy. It is the dominant part when the interaction |cu| of the 
states is small. The value lm{Z) is dominant when |u;| is large. It is related, according to 
( 1A2I) . to a bifurcation of the widths of the levels. 

The formal equivalence of the optical wave equation in PT symmetric optical lattices to 
the quantum mechanical Schrodinger equation allows us to study the properties of quantum 
systems the states of which can not only decay due to their cou pling to the environment 
according to (1A3I) . but may also be formed out of the environment [2J In optics, these two 
oossibilities are called loss and gain. In PT symmetric optical lattices, the eigenvalues are 

121 

£1,2 = ei, 2 ± 2 7i,2 ( wit h 7i,2 > and e 1 = e 2 ) (A4) 
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in difference to ( 1A3j) . Due to PT symmetry, all eigenvalues e, = e« may be real (correspond- 
ing to 7; = 0) when Re(Z) 3> Im(Z), i.e. at low coupling of the states to the continuum. 
However, the PT symmetry is broken when Im(Z) 3> Re(Z) and 7^2 7^ 0. The difference 
between the two models with the eigenvalues fl A3j) and (IA4[) will allow us to receive inter- 
esting information on quantum systems by studying not only open quantum systems (which 
exist in nature) but also PT symmetric systems (which are formally equivalent to them). 



Appendix B: The eigenfunctions of a non-Hermit ian 2x2 operator 

The eigenfunctions of the non-Hermitian Hamilton operator H, equation (I A 1 1) . are 
biorthogonal, 

(<j>*M = Sk,i ■ (Bl) 

From these equations follows 



{<Pk\<Pk) = A k > 1 (B2) 



B l k ; \B l k \>0. (B3) 



At the crossing point A k — > 00, \B l k ^| — > 00, for details see 

The relation between the eigenfunctions 0i and 02 of the operator ( 1A1I) at the crossing 
point is 

07 ± % ^ ; c 2 r T * 0f (B4) 

according to analytical as well as numerical studies That means, the state 0i jumps, at 
the exceptional point, via the chiral state <pi ± i<p2 to the state ±i0 2 - 

The two eigenfunctions are linearly dependent of one another at the crossing point such 
that the number of eigenfunctions of H seems to be reduced at this point. Theoretical 
studies (4?]] have shown however that associated vectors 4>^ ra defined by the Jordan relations, 
appear at the crossing points. The corresponding equations are 

(H - £ ) C 2 = ; (H - s ) ft™ = 0^ • (B5) 

The existence of two states in the very neighborhood of the exceptional point has been seen 



in a numerical calculation for the elastic scattering of a proton on a light nucleus 
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48|: the 



elastic scattering phase shifts jump always by 2tt (and not by tt as for a single resonance 
state). 

n 

In an experimental study on a microwave cavity [46|, the topological structure of the 
exceptional point and its surrounding is studied by encircling it and tracing the relative 
amplitudes of the wavefunctions (field distributions inside the cavity). As a result, the 
wavefunctions including their phases are restored after four surroundings. The authors 46] 
interpreted the experimental data by two theoretical assumptions: (i) the two wavefunctions 
coalesce into one at the exceptional point, 0f -H- 0f, and (ii) only one of the wavefunctions 
picks up a phase of tt (a sign change) when encircling the critical point. The experimental 
result can be explained, however, without any additional assumptions by using the relations 
(ED: 

1. cycle: e 1)2 ->■ e 2 ,i; 0i, 2 ± i <f>2,i ] 2. cycle: e 2 ,i ->• e lj2 ; ± % 2 ,i ->> -0i, 2 
3. cycle: e 1>2 e 2) i; -0 1)2 ->■ ; 4. cycle: e 2)1 ->■ e 1>2 ; =Fi0 2) i->- 0i, 2 . 

The eigenvalues are restored after two surroundings and the eigenfunctions are restored after 









1 <I>1 1 







at the crossing point. The topological phase is twice the Berry phase, in accordance with 
the enlarged function space in open quantum systems. 

Furthermore, the phases of the wavefunctions jump by 7r/4 at the crossing point (when 
traced as a function of a parameter) due to the biorthogonality fIBlj) of the eigenfunctions 
of the non-Hermitian Hamiltonian H, see also ( 1B2I) and (IB 3 1) . This result has been proven 
in many numerical studies, see {§(]. 



Appendix C: The phase rigidity r\ of the eigenfunctions of a non-Hermitian 2x2 
operator 

Let us now consider the consequences of the biorthogonality relations fIBll) and (IB2I) for 
the two borderline cases characteristic of neighboring resonance states. 

(i) The two levels are distant from one another. Then the eigenfunctions are (almost) 
orthogonal, (<Pl\4> k ) ~ {4>k\4>k) = A h w 1. 

(ii) The two levels cross. Then the two eigenfunctions are linearly dependent according to 
dE3j) and {(j) k \(j) k ) =A k ^oo. 

These two relations show that the phases of the two eigenfunctions relative to one another 
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change when the crossing point is approached. This can be expressed quantitatively by 
defining the phase rigidity of the eigenfunctions <pk, 

r. - = V ■ (CI) 

It holds 1 > r/j > 0. The non-rigidity of the phases of the eigenfunctions of H follows 
also from the fact that (4>l\4>k) is a complex number (in difference to the norm (4>k\4>k) which 
is a real number) such that the normalization condition ( IB 11) can be fulfilled only by the 
additional postulation lm{(f>* k \(p k ) = (what corresponds to a rotation). 

Note, discrete states never cross but always avoid crossing, see appendix [A] The eigen- 
functions are real and normalized as (4>k\4>k) — 1 (see appendix |A|) . Accordingly, r k = 1 for 
discrete states also in the region of an avoided level crossing. The corresponding crossing 
(exceptional) point can be found by analytical continuation into the continuum [181 ] . 

The variation of r k in approaching the crossing point of two eigenvalue trajectories of 
resonance states is proven experimentally by means of a study on a microwave cavity 49]. 
As a result, the phase difference between two modes is tc at large distance and decreases to 
7r/2 at the crossing point. The authors of [49j interpret the experimental data by assuming 
(i) that the singular point is a chiral state (in spite of the phase jump occurring at the 
crossing point, when traced as a function of a certain parameter, see ( 1B4I) ). (ii) that the 
number of states is reduced from 2 to 1 at the crossing point (in spite of the existence of 
the associate vector flB5j) ) and (iii) that a single point in the continuum can be identified 
(although it is of measure zero). The authors are unable to explain the large parameter 
range in which the phase difference decreases in approaching the crossing point. 

Considering the phase rigidity r k in the regime of the two overlapping resonance states, 
no additional assumptions are required for the explanation of the experimental results given 
in 49J], since the phase rigidity (being a quantitative measure for the degree of resonance 
overlapping) varies smoothly in a comparably large parameter range. It can therefore be 
concluded that the experimental results prove the statement that the phases of the 
eigenfunctions of the non-Hermitian Hamilton operator H, equation (1A1I) . are not rigid in 
approaching the crossing point. 
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Appendix D: Nonlinear source term in the Schrodinger equation in the neighbor- 
hood of an exceptional point 



According to (IA1B . the Schrodinger equation with the unperturbed operator Hp = H(cj 
0) and a source term arising from the interaction u with another state reads 50] 



(H - e n ) \<j> n ) = - \d> n ) = W \<f> n ) 




= E (<Pk\W\<Pn){A k \4> k ) +J2 B l k . (Dl) 

fc=l,2 l^k 

Here (<p k \<Pk) = A k > 1 according to (|B2]) and (0 fc |0 ¥A .) = -(fa^fa) = B l k , \B l k \ > 
according to (1B3I) . The A k and B l k characterize the degree of resonance overlapping. In the 
regime of overlapping resonances, 1 > A k > 0, \B k \ > 0, and equation (IDlj) is nonlinear. 
The most important part of the nonlinear contributions is contained in 

(H - e n ) \<j) n ) = {<l> n \W\<t> n ) |0„| 2 \(j> n ) (D2) 

which is a nonlinear Schrodinger equation. According to (1D1I) . the nonlinear Schrodinger 
equation (ID2I) passes smoothly into the standard linear Schrodinger equation when A k — > 1 
and B{ 0. 



Appendix E: Time reversal symmetry breaking in the neighborhood of an excep- 
tional point 

Exceptional points that are well separated from the influence of external sources (includ- 
ing the influence caused by other resonance states), are highly symmetric in approaching 
them. That means, the two states pass one into the other one according to (1B4|) with an 
exchange of their wavef unctions, </>f — > ± i <pf and 0f — > =F i 0f. At a certain finite 
distance from the exceptional point, there are again two states with the wavefunctions \<f>i\ 
and 1 4>2 1, respectively. 

This symmetry may be distorted under the influence of an external magnetic field as 



51| . The magnetic field causes time 



has been shown experimentally on a microwave cavity 
reversal symmetry breaking. 

The symmetry may be disturbed also by the influence of another resonance state in the 
neighborhood due to the finite parameter range around the exceptional point in which the 
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501 ] . When the interaction of the 



wavefunctions of the two states are mixed with each other 
third state is symmetric relative to the two crossing ones, the third state will appear as an 
observer and time reversal symmetry is not broken. Numerical examples of such a situation 
arc shown in the transmission through a quantum dot 
in 53] . When the interaction of the third state with the two crossing ones is, however, not 



52] and also in the generic case studied 



symmetrically, time reversal symmetry may be broken and may cause irreversible processes 
due to the nonlinear terms in the Schrodinger equation as discussed in appendix [Dl 

It would be highly interesting to study experimentally time reversal symmetry breaking 
in the case of an exceptional point disturbed non-symmetrically by a third state in the 
neighborhood. 
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